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The results of a comparative analysis based upon a Karhunen-Loeve expansion of turbulent pipe 
flow and drag reduced turbulent pipe flow by spanwise wall oscillation are presented. The turbulent 
flow is generated by a direct numerical simulation at a Reynolds number Re T = 150. The spanwise 
wall oscillation is imposed as a velocity boundary condition with an amplitude of A + — 20 and a 
period of T + = 50. The wall oscillation results in a 27% mean velocity increase when the flow is 
driven by a constant pressure gradient. The peaks of the Reynolds stress and root-mean-squared 
velocities shift away from the wall and the Karhunen-Loeve dimension of the turbulent attractor is 
reduced from 2453 to 102. The coherent vorticity structures are pushed away from the wall into 
higher speed flow, causing an increase of their advection speed of 34% as determined by a normal 
speed locus. This increase in advection speed gives the propagating waves less time to interact 
with the roll modes. This leads to less energy transfer and a shorter lifespan of the propagating 
structures, and thus less Reynolds stress production which results in drag reduction. 



I. INTRODUCTION 

In the last decade a significant amount of work 
has been performed investigating the structure of wall 
bounded turbulence, with aims of understanding its self- 
sustaining nature and discovering methods of control. 1 
One of the greatest potential benefits for controlling tur- 
bulence is drag reduction. As the mechanics of the 
different types of drag reduction are studied, most ex- 
planations of the mechanism revolve around control- 
ling the streamwise vortices and low speed streaks j 2 -* 3 - 
One such method of achieving drag reduction is span- 
wise wall oscillation, first discovered by Jung et ali^ in 
1992, and later confirmed both numerically^^^ and 
experimentally ^ 10 ! 11 ! 12 ! 13 to reduce drag on the order of 
45%. The prevalent theory of the mechanism behind this 
was developed using direct numerical simulation of a tur- 
bulent channel flow by Choi et al. 6 and experimentally 
confirmed by Choi and Clayton^ showing that the spa- 
tial correlation between the streamwise vortices and the 
low speed streaks are modified so that high speed fluid 
is ejected from the wall, and low speed fluid is swept to- 
wards the wall. Even though this proposed mechanism 
describes the near- wall dynamics that govern the drag re- 
duction, questions behind the global dynamics have not 
been sufficiently resolved. What is the effect (if any) on 
the outer region of the flow? How do the coherent struc- 
tures of the flow near the wall adjust with the spanwise 
oscillations? What is the effect on the interactions be- 
tween the inner and outer layers? 

One manner in which to address these questions is 
through direct numerical simulation (DNS) of turbu- 
lence. As supercomputing resources increase, DNS con- 
tinues to provide an information rich testbed to inves- 
tigate the dynamics and mechanisms behind turbulence 
and turbulent drag reduction. DNS resolves all the scales 
of turbulence without the need of a turbulent model and 
provides a three dimensional time history of the entire 
flow field. One of the methods used for mining the 



information generated by DNS is the Karhunen-Loeve 
(KL) decomposition, which extracts coherent structures 
from the eigenfunctions of the two-point spatial corre- 
lation tensor. This allows a non-conditionally based in- 
vestigation that takes advantage of the richness of DNS. 
The utility of this method is evident in the knowledge 
it has produced so far, such as the discovery of propa- 
gating structures (traveling waves) of constant phase ve- 
locity that trigger bursting and sweeping events i 15 ! 16 ! 17 
These studies in turn have lead to a new class of meth- 
ods for achieving drag reduction through wall imposed 
traveling waves Another study using KL decomposi- 
tion examined the energy transfer path from the applied 
pressure gradient to the flow through triad interaction of 
structures j 19 i 2Q explaining the dynamical interaction be- 
tween the KL modes. In the realm of control, KL meth- 
ods have been used to produce drag reduction in a turbu- 
lent channel by phase randomization of the structures 2 ^ 
and to understand the effect of drag reduction by con- 
trolled wall normal suction and blowing. 22 In the present 
study, the KL framework is used to examine the differ- 
ences in the turbulent structures and dynamics between 
turbulent pipe flow with and without spanwise wall os- 
cillation. 

For this comparative analysis, turbulent pipe flow was 
chosen as opposed to turbulent channel flow because of its 
industrial relevance and experimental accessibility. The 
main difference between pipe and channel flows is that in 
turbulent pipe flow the mean flow profile exhibits a log- 
arithmic profile that overshoots the theoretical profile at 
low Reynolds numbers, whereas in turbulent channel flow 
it does not i 23 i 24 Secondly, pipe flow differs from chan- 
nel flow because pipe flow is linearly stable to infinitesi- 
mal disturbances where channel flow is linearly unstable 
above a critical Reynolds number j 25 i 26 A significant com- 
putational difference between pipe and channel flow is 
the presence of a numerical difficulty introduced by the 
singularity in polar-cylindrical coordinates at the pipe 
centerline, which has limited the number of DNS studies 
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in turbulent pipe flow i 27 ' 28 i 29 i 3Q ' 31 i 32 Similarly, many KL 
studies have been performed in a turbulent channel flow, 
but to the best of our knowledge the work presented here 
is the first to extend the KL method to spanwise wall os- 
cillated turbulent pipe flow. 

In previous work a DNS of turbulent pipe flow 
for Re T = 150 was benchmarked and its KL ex- 
pansion was reported, forming the baseline for this 
study. 33 Similar structures to those of turbulent channel 
flo w 15 i 16 i 17 were found, including the presence of propa- 
gating modes. These propagating modes are character- 
ized by a nearly constant phase speed and are responsible 
for the Reynolds stress production as they interact and 
draw energy from the roll modes (streamwise vortices). 15 
Without this interaction and subsequent energy trans- 
fer, the propagating waves decay quickly, reducing the 
total Reynolds stress of the flow #22 As shown by Sirovich 
et al. ^ 15 i 16 the interaction between the propagating waves 
and the roll modes occurs by the propagating waves form- 
ing a coherent oblique plane wave packet. This wave 
packet interacts with the roll modes, and when given 
enough interaction time, the roll mode is destabilized 
eventually resulting in a bursting event. 15 It is in this 
bursting event that the energy is transferred from the 
rolls to the propagating waves. 19 In this paper we show 
that in the presence of spanwise wall oscillation these 
propagating modes are pushed away from the wall into 
higher speed flow. This causes the propagating modes 
to advect faster, giving them less time to interact with 
the roll modes. This leads to reduced energy transfer 
that occurs less often, and yields lower Reynolds stress 
production, which ultimately results in drag reduction. 

II. NUMERICAL METHODS 

We use a globally high order spectral element Navier- 
Stokes algorithm to generate turbulent data for pipe flow 
driven by a mean streamwise pressure gradient i 35 i 36 The 
non-dimensional equations governing the fluid are 

d t U + U • VU = - VP + Re" 1 V 2 U (1) 
V-U = 0, (2) 

where U is the velocity vector, Re r is the Reynolds 
number, and P is the pressure. The velocity is non- 
dimensionalized by the wall shear velocity U T = \/r w /p 
where r w is the wall shear stress and p is the density. 
The Reynolds number is Re r = U T R/v = 150, where R 
is the radius of the pipe, and v is the kinematic viscosity. 
When non-dimensionalized with the centerline velocity, 
the Reynolds number is Re c w 4300. Two simulations 
were performed, one with and one without spanwise wall 
oscillation. In a pipe, the spanwise direction corresponds 
to the azimuthal direction, so the oscillation is about the 
axis of the pipe. Each case was run for t + — U^t/v = 
16800 viscous time units. In the oscillated case, the sim- 
ulation was performed with an azimuthal velocity wall 




FIG. 1: A cross-section of a spectral element grid used for the 
pipe flow simulation. The arrow denotes the direction of the 
spanwise (azimuthal) oscillation about the axis of the pipe. 



boundary condition vq(t = R, 9, z) = A + sin(27rt/T + ) of 
amplitude A+ = A/U T = 20 and period T+ = U*T/v = 
50 with (r, z) being the radial, azimuthal, and stream- 
wise coordinates respectively. This is not intended to be 
a parametric study and the amplitude and period were 
chosen to achieve maximum possible drag reduction be- 
fore relaminarization occurred. 

To solve equations [T] and [2] we use a numerical al- 
gorithm employing a geometrically flexible yet expo- 
nentially convergent spectral element discretization in 
space. The spatial domain is subdivided into elements, 
each containing a high-order (12th order) Legendre La- 
grangian interpolanti^ 7 . The spectral element algorithm 
elegantly avoids the numerical singularity found in polar- 
cylindrical coordinates at the origin, as seen in Figure 
[TJ The streamwise direction contains 40 spectral ele- 
ments over a length of 10 diameters. The effective res- 
olution of the flow near the wall is Ar + w 0.78 and 
(RA6)+ w 4.9, where the radius r + and the arc length at 
the wall (RA6)+ are normalized by wall units v/U T de- 
noted by the superscript +. Near the center of the pipe, 
the grid width is A + w 3.1. The grid spacing in the 
streamwise direction is a constant Az + = 6.25 through- 
out the domain. Further details can be found in Duggleby 
et ali22 

The flow is driven by a constant mean pressure gradi- 
ent to keep Re r constant. The spanwise wall oscillation 
results in a mean flow rate increase, effectively changing 
the Reynolds number based upon mean velocity (Re m ) 
while keeping Re r constant. This keeps the dominant 
structures of the flow similar, as they are affected primar- 
ily by the inner layer wall shear stress. 7 The oscillations 
were started on a fully turbulent pipe at Re r = 150, and 
to avoid transient effects, data were not taken until the 
mean flow rate had settled at its new average value over 
a time interval of 1000t + . 

In the Karhunen-Loeve (KL) procedure, the eigenfunc- 
tions of the two-point velocity correlation tensor, defined 

by 
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/ / / K(x,x')<f>(x')r'dr'de'dz' = (3) 
Jo Jo Jo 

K(x,x / ) = (u(x)0u(x / )), (4) 

are obtained, where x = (r, 0, z) is the position vector, 
3>(x) is the eigenfunction with associated eigenvalue A, 
if (x, x') is the kernel, and (g) denotes an outer product. 
In order to focus on the turbulent structures, the kernel 
is built using fluctuating velocities u = U — U. The 
mean velocity, U, is found by averaging over all #, z, and 
time. The angle brackets in equation [4] represent the time 
average using an evenly spaced time interval over a total 
time period sufficient to sample the turbulent attractor. 
In this study, the flow field was sampled every 8t + for a 
total time of 16800£+. 

Since the azimuthal and streamwise directions are pe- 
riodic, the kernel in the azimuthal and streamwise direc- 
tion is only a function of the distance between x and x' 
in those respective directions. Therefore the kernel can 
be rewritten as 

K(r,6,z,r',e',z') = K(r,r \0 - 0' ', z - z) 

= JC(m,n',r,r f )e in9 e i27Tmz/L (5) 

with azimuthal and streamwise wavenumbers n and m 
respectively and the remaining two-point correlation in 
the radial direction /C(m, n; r, r'). It can be shown that 
in this form, the Fourier series is the resulting KL func- 
tion in the streamwise and azimuthal direction. 38 The 
resulting eigenfunction then takes the form 

*(r, 0, z) = *(m, n; r)e m0 e^ /L . (6) 

Making use of this result, and noting that the two-point 
correlation in a periodic direction is simply the Fourier 
transform of the velocities, the azimuthal and stream- 
wise contributions to the eigenfunctions are extracted a 
priori by taking the Fourier transform of the velocities 
u(r, 0, z) = J2m n=o u(ra, n \ r)e m6> e z27rm2 / L and forming 
the remaining kernel /C(m,n;r, r') for each wavenumber 
pair n and m. The eigenfunction problem, with the or- 
thogonality of the Fourier series taken into account, is 




/C(m, n; r, r')\l>* (r')r' ' dr' = A mn \l>(ra, n; r), (7) 



/C(m, n; r, r') = (u(m, n; r) u*(ra, n; r')). (8) 

where the ★ denotes the complex conjugate since the func- 
tion is now complex, and the weighting function r' is 
present because the inner product is evaluated in polar- 
cylindrical coordinates. The final form is still Hermitian 
just as it was in equation [3j the discrete form of equation 
[3 is kept Hermitian by splitting the integrating weight 
and solving the related eigenvalue problem 

[y/r^lCp 8 (m, n; r p , r s ) y/rl] [^/r~ s ^ q {m, n; r a )]* 

= ^mnq [^^f q {m, n, ; r p j\ , (9) 



where /C ps (m, n; r p , r 3 ) is the discretization of 
/C(m, n; r, r') using a Q point quadrature to evalu- 
ate equation [71 with p, s = 1,2, Because the 
kernel is built with the two-point correlation between all 
three coordinate velocities, its solution has 3Q complex 
eigenfunctions *S? q and corresponding eigenvalues, listed 
in decreasing order \ mnq > ^ m n(q+i) f° r a given m and 
n, with quantum number q = 1,2, ...,3Q. It is noted 
that equation [9] is only valid for a trapezoidal integration 
scheme with evenly spaced grid points (which was used 
in this study), a different quadrature with weight w(r) 
can be incorporated in a similar fashion, keeping the 
final matrix Hermitian. 

The eigenfunctions\l/ g (m, n; r) hold certain properties. 
Firstly, they are normalized with inner product of unit 
length tyq^q/rdr = 5 qq /, where S is the Kronecker 
delta. Secondly, since the eigenfunctions represent a flow 
field, 

^q(m,n;r) = (\I/^(m, n; r), ^g(m, n; r), ^g(m, n; r)) T 

(10) 

with radial, azimuthal, and streamwise components 
^^(m,n;r), \I/^(ra,n;r), and ^^(m,n;r) respectively, 
they hold the properties of the flowfield such as boundary 
conditions (no slip) and continuity 

1 d . _ _ , >. in t9 , N i2n:m _ z/ , 

~ — (rw q (m, n; r))-\ w q [m, n; r)H — n; r) = 0. 

r dr r 1j 

(ii) 

Thirdly, the eigenvalues represent the average energy of 
the flow contained in the eigenfunction ^(m, n; r) 

Xmn q = ( | u(r, 9, z) , * q (m, ft; r)e ine e i2 ™ z ' L \ 2 ) , (12) 

which is why it is necessary that the discrete matrix in 
equation [9] must be Hermitian to get real and positive 
eigenvalues. These three properties allow the eigenfunc- 
tions to be tested for their validity. 

In summary, the KL procedure yields an orthogonal 
set of basis functions (modes) that are the most energet- 
ically efficient expansion of the flow field. By studying 
the subset that includes the largest energy modes, in- 
sight is gained as this subset forms a low dimensional 
model of the given flow. Examining the structure, dy- 
namics, and interactions of this low dimensional model 
yields important information that we use to build a bet- 
ter understanding of the dynamics of the entire system. 



III. RESULTS 

Spanwise wall oscillation results in four main effects on 
the flow, its structures, and its dynamics. They are: 

1. An increase in flow rate and a shifting away from 
the wall of the root-mean-squared velocities and 
Reynolds stress peaks. 
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FIG. 2: Mean velocity fluctuations (solid) and the average 
velocity (dash-dot) for the oscillated case versus time (£ + ). 
The fluctuations are consistent with turbulent flow, and the 
mean velocity is 26.9% greater than in the non-oscillated pipe. 



2. A reduction in the dimension of the chaotic attrac- 
tor describing the turbulence. 

3. An increase in energy of the propagating modes re- 
sponsible for carrying energy away from the wall to 
the upper region, while the rest of the propagating 
modes exhibit a decrease in energy. 

4. An increase of the advection speed of the traveling 
wave packet as determined by a normal speed locus. 

First, surface oscillation has a major effect on the tur- 
bulent statistics. The spanwise wall oscillation of ampli- 
tude A + = 20 and period T + = 50 resulted in a flow 
rate increase of 26.9%, shown in Figure [2j This com- 
bination of amplitude and period was chosen because it 
provides the largest amount of drag reduction while keep- 
ing the flow turbulent. Numerical simulations with larger 
oscillations completely relaminarize the flow. The com- 
parison of the mean profiles in Figures [3] and 2] show the 
higher velocity in the outer region, with the inner region 
remaining the same, as expected by keeping a constant 
mean pressure gradient across the pipe. 

The comparison of the root-mean-square (rms) veloc- 
ity fluctuation profiles and the Reynolds stress profile in 
Figures [5] and [6] show that the streamwise fluctuations 
decrease in intensity by 7.5% from 2.68 to 2.48. Also, 
the change in peak location from y + = 16 to y + = 22 
away from the wall has the same trend as the maxi- 
mum Reynolds stress u r u z , where y + = (R — r)U T / v is 
the distance from the wall using normalized wall units 
(y/U T ). The azimuthal fluctuating velocities show a 
slightly greater magnitude peak of 1.03 closer to the wall 
at y + = 29 versus 0.99 at y + = 40 for the non-oscillated 
pipe. The radial fluctuations remain almost unchanged, 
showing only a slight decrease from the wall through the 
log layer (y + « 100), resulting in the peak shifting from 
0.81 at y+ = 55 to 0.78 at y+ = 61. The Reynolds stress 



FIG. 3: Mean velocity profile for a non-oscillated (solid) and 
oscillated (dashed) turbulent pipe flow versus y + . Theoretical 
(dash-dot) includes the sublayer (u + — y + ) and the log layer 
(u + = log(?/ + )/0.41+5.5). The mean profile shows a log layer, 
but overshoots the theoretical value as expected for pipe flow 
until a much higher Reynolds number. 
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FIG. 4: Mean velocity profile for non-oscillated (solid) and 
oscillated (dashed) pipes with their respective bulk velocities 
(dash- dot and dots) versus radius show an increase in bulk 
velocity of 26.9%. 



also shows a reduction in strength and a shift away from 
the wall. The peak changes from 0.68 at y + = 31 to 0.63 
at y + = 38. Thus, a major difference between the two 
flow cases, in addition to the expected flow rate increase, 
is the shift of the rms velocity and Reynolds stress peaks 
away from the wall. 

The second major effect can be found by examining 
the size of the chaotic attractor describing the turbu- 
lence. The eigenvalues of the KL decomposition repre- 
sent the energy of each eigenfunction. By ordering the 
eigenvalues from largest to smallest, the number of eigen- 
functions needed to capture a given percentage of energy 
of the flow is minimized. Table [J shows the 25 most en- 
ergetic eigenfunctions, and Figure [71 shows the running 
total of energy versus mode number. 90% of the energy 
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FIG. 5: Root-mean-squared velocity fluctuations for the non- 
oscillated (solid) and oscillated (dashed) case versus y + . The 
oscillated case shows a shift away from the wall, except for 
the azimuthal rms which captures the Stokes layer imposed 
by the wall oscillation. 



FIG. 7: Comparison of the running total energy retained in 
the KL expansion for the non-oscillated (solid) and oscillated 
(dashed) cases. The 90% crossover point is 2453 and 102 
respectively. This shows a drastic reduction in the dimension 
of the chaotic attractor. 
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FIG. 6: Reynolds Stress u r u z versus y + for non-oscillated 
(solid) and oscillated (dashed) cases. Similar to the rms ve- 
locities, the Reynolds stress shows a shift away from the wall 
from y + — 31 to y + — 38. 



FIG. 8a: (0,0,1) eigenfunction capturing the imposed velocity 
oscillations. The radial velocity is zero except near the wall, 
where a slight non-zero component is evident, induced by the 
spanwise oscillation. 



is reached with Dkl = 102 compared to Dkl = 2453 
for the non-oscillated case. This mark, known as the 
Karhunen-Loeve dimension, is a measure of the intrin- 
sic dimension of the chaotic attractor of turbulence as 
discussed by Sirovic h 39 i 4Q . By oscillating the wall our re- 
sults show that the size of the attractor is reduced, and 
the system is less chaotic. 

The third major effect is found by examining the en- 
ergy of the eigenfunctions. The 25 most energetic are 
listed for each case in Table HI The top ten modes 
with the largest change in energy are shown in Table HI1 
Firstly, the order of the eigenfunctions remain relatively 
unchanged, with a few notable differences. The (0,0,1) 
and (0,0,3) shear modes represent the Stokes flow as seen 
in Figures [8a| and [8cl and the (0,0,2) mode, shown in 
Figure l8b( represents the changing of the mean flow rate, 



similar to the non-oscillated (0,0,1) mode. The (1,2,1) 
mode shows a large increase in energy. Also of note is 
the reduction in energy of the (3,1,1) and (4,1,1) modes, 
but their structure remains virtually unchanged. 

In examining the energy content of the structure sub- 
classes as a whole, a trend is discovered, shown in Ta- 
ble IIIII Each of these subclasses, reported in Duggleby 
et al., 33 were found to have similar qualitative coher- 
ent vorticity structure associated with their streamwise 
and azimuthal wave number. We use "coherent vortic- 
ity" to refer to the imaginary part of the eigenvalues 
of the strain rate tensor duijdx^ following the work of 
Chong et al.. 41 Based upon the qualitative structure, the 
propagating or traveling waves, described by non-zero 
azimuthal wavenumber and nearly constant phase speed, 
were found to have four subclasses: the wall, the lift, the 
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TABLE I: The 25 most energetic 
wavenumber, n is the spanwise 
eigenvalue quantum number. 



modes; m is the streamwise 
wavenumber, and q is the 



Non-oscillated 
Index m n q Energy % Total 



Oscillated 
m n q Energy % Total 
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0.48% 


8 


1 


3 


1 
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FIG. 8b: (0,0,2) eigenfunction, similar in structure to the 
non-oscillated (0,0,1). 
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FIG. 8c: (0,0,3) eigenfunction capturing the secondary flow 
in the Stokes layer. This imposed Stokes flow created by the 
oscillation dominates the near wall region of the oscillated 
pipe. 



TABLE II: Ranking of eigenfunctions by energy change be- 
tween the non-oscillated and the oscillated cases, m is the 
streamwise wavenumber, n is the spanwise wavenumber, and 
q is the eigenvalue quantum number. 



Increase Decrease 
Rank AAk m n q AAk m n q 
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215 








1 


-0.962 
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34.5 
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1 
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asymmetric, and the ring modes. 

The wall modes are found when the spanwise 
wavenumber is larger than the streamwise wavenum- 
ber. They possess a qualitative structure having co- 
herent vortex cores near the wall, and have their en- 
ergy decreased by 20.4% with wall oscillation. Likewise, 
the ring modes, which are found for non-zero streamwise 
wavenumber and zero spanwise wavenumber with rings 
of coherent vorticity, have their energy decreased by 5.2% 
with wall oscillation. The asymmetric modes have non- 
zero streamwise wavenumber and spanwise wavenumber 
n = 1, which allows them to break azimuthal symmetry. 
These modes undergo a decrease in energy of 2.3% with 
spanwise wall oscillation. 

Conversely to the decrease in energy found in the other 
three propagating modes, the lift modes increase in en- 
ergy by 1.5% with spanwise wall oscillation. These modes 
are found with a streamwise wavenumber that is greater 
than the spanwise wavenumber, and they display coher- 
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TABLE III: Energy comparison of turbulent pipe flow struc- 
ture subclasses between non-oscillated and oscillated pipes, m 
is the streamwise wavenumber, and n is the azimuthal (span- 
wise) wavenumber. All the propagating modes decrease in 
energy, except the lift modes. 



Structure 


Energy 

Non-oscillated Oscillated 


Propagating Modes (m > 0) 


53.48 


48.86 


(a) Wall (n > m) 


23.5 


18.7 


(b) Lift (m > n) 


19.8 


20.1 


(c) Asymmetric (n — 1) 


6.07 


5.93 


(d) Ring (n = 0) 


4.41 


4.18 


Non-propagating Modes (m = 


0) 12.97 


265 


(a) Roll mode (n > 0) 


12.3 


13.2 


(b) Shear mode (n = 0) 


0.72 


251.8 



ent vortex structures that starts near the wall and lifts 
away from the the wall to the upper region. Combined, 
the four propagating subclasses modes lose 9.02% of their 
energy, whereas the non-propagating modes (the modes 
with zero azimuthal wavenumber) gain 2043%. 

Thus, the third major result is that the energy of the 
propagating wall, ring, and asymmetric modes decrease 
while the energy of lift mode increases slightly. Fol- 
lowing the work of Sirovich et al^ this shows that en- 
ergy transfer from the streamwise rolls to the traveling 
waves is reduced, and any energy that is transferred is 
quickly moved away from the wall to the outer region 
(lift modes). The energy spectra showing the change of 
energy by subclasses is shown in Figure [U 




Streamwise Wavenumber 



FIG. 9: Comparison of energy spectra for non-oscillated 
(solid) and oscillated (dashed) flows for the propagating mode 
subclasses. All propagating subclasses in the oscillated case 
decrease with respect to the non-oscillated case, except for 
the lift mode, which increases slightly. 



The fourth and most important effect is that the prop- 
agating modes advect faster in the oscillated case. The 
normal speed locus of the 50 most energetic modes of 
both cases is shown in Figure [10l For this, the phase 
speed cj/||k|| is plotted in the direction k/||k|| with k = 
(m, n). A circular locus is evidence that these structures 
propagate as a wave packet or envelope that travels with 
a constant advection speed. The advection speed is given 
by the intersection of the circle with the abscissa. By ex- 
amining the normal speed locus of the non-oscillated and 
oscillated pipe flow, the wave packet shows an increase in 
advection speed from 8.41 U T to 10.96Z7 r , an increase of 
30%. This is a result of the oscillating Stokes layer push- 
ing the structures away from the wall into a faster mean 
flow by creating a dominant near wall Stokes layer where 
the turbulent structures cannot form. This is confirmed 
by the shifting of the rms velocities and Reynolds stresses 
away from the wall as reported earlier. In addition to 
a faster advection speed, the energy of the propagating 
modes decay faster, resulting in bursting events with a 
shorter lifespan. This is seen in Figure [TT1 where the av- 
erage burst duration of the (1,5,1) mode is reduced from 
106t + in the non-oscillated case to 65.3£ + in the oscil- 
lated case. The burst duration is taken to be the average 
time of all events where the square of the amplitude of 
the mode is more than one standard deviation greater 
than the mean. 
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FIG. 10: Comparison of the normal speed locus for the oscil- 
lated (•) and non-oscillated (+) case. The solid and dashed 
lines represent a circle of diameter 8.41 and 10.96 respectively 
that intersect at the origin. 

The shifting of the structures away from the wall is 
shown for the most energetic modes for each propagat- 
ing subclass in Figures [T2l through [2T1 These structures, 
which represent the coherent vorticity of the four sub- 
classes of the propagating modes, are pushed towards the 
center of the pipe, where the mean flow velocity is faster. 
The location of the coherent vortex core for these eight 
modes are listed in Table [IVl showing a shift away from 
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FIG. 11: A reduction in the burst duration of the (1,5,1) 
mode from 106t + for the non-oscillated case (top) to 65.3£ + 
in the oscillated case (bottom) shows a faster decay of the 
bursting energy with spanwise wall oscillation. The burst 
duration is the average time of all events where the square of 
the amplitude (|a(£)| 2 ) is more that one standard deviation 
greater than the mean. This amplitude level is denoted by 
the dashed line. 
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FIG. 12: Cross-section of coherent vorticity of (1,2,1) wall 
mode. Left: (a) non-oscillated. Right: (b) oscillated. The 
vortex core shifts y + = 6.8 away from the wall. 
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FIG. 13: Cross-section of coherent vorticity of (1,5,1) wall 
mode. Left: (a) non-oscillated. Right: (b) oscillated. The 
vortex core shifts y + — 9.9 away from the wall. 
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FIG. 14: Cross-section of coherent vorticity of (2,2,1) lift 
mode. Left: (a) non-oscillated. Right: (b) oscillated. The 
vortex core shifts y + — 11.0 away from the wall. 



and that they had higher phase velocities. Thirdly, the 
study by Zhou 7 is also consistent with these results, as 
she found that any oscillation in the streamwise direc- 
tion reduces the effectiveness of the drag reduction. Any 
streamwise oscillation, even though it would still push 
the structures away from the wall, would adversely effect 
the mean flow rate profile resulting in an advection speed 
of the propagating waves that is less than in the purely 
spanwise oscillated case. 



the wall. The only mode found not to follow this trend is 
the (1,0,1) mode, which undergoes a major restructuring 
resulting in its vortex core moving towards the wall. 

This faster advection explains the experimental results 
found by K.-S. Choi 42 that showed a reduction in the 
duration and strength of sweep events in a spanwise wall 
oscillated boundary layer of 78% and 64% respectively. 
For this experiment, the flowrate was kept constant, so 
the energy was reduced, whereas in our case the mean 
pressure gradient was kept constant yielding virtually no 
change in the energy of the propagating structures. Also 
corroborating these results is the work by Prabhu et al*22 
that examined the KL decomposition of controlled suc- 
tion and blowing to reduce drag in a channel. They, too, 
found that the structures were pushed away from the wall 
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FIG. 15: Cross-section of coherent vorticity of (3,2,1) lift 
mode. Left: (a) non-oscillated. Right: (b) oscillated. The 
vortex core shifts y + — 11.2 away from the wall. 
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FIG. 16: Cross-section of coherent vorticity of (1,1,1) asym- 
metric mode. Left: (a) non-oscillated. Right: (b) oscillated. 
The vortex core shifts y + — 3.2 away from the wall. 
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FIG. 17: Cross-section of coherent vorticity of (2,1,1) asym- 
metric mode. Left: (a) non-oscillated. Right: (b) oscillated. 
The vortex core shifts y + — 45.9 away from the wall with a 
slight change in structure. 
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FIG. 20: Cross-section of coherent vorticity of (0,6,1) roll 
mode. Left: (a) non-oscillated. Right: (b) oscillated. The 
vortex core shifts y + — 10.1 away from the wall. 





Qwerem VTOfliCHyi qjQE Q.0B 0.1 0.1? 





0.03 0J04 0.05 0.05 0J07 



FIG. 21: Cross-section of coherent vorticity of (0,2,1) roll 
mode. Left: (a) non-oscillated. Right: (b) oscillated. The 
vortex core shifts y + = 7.8 away from the wall. 
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FIG. 18: Cross-section along the r-z plane of coherent vor- 
ticity of (1,0,1) ring mode. Top: (a) non-oscillated. Bottom: 
(b) oscillated. The vortex core shifts y + — 13 towards the 
wall with significant changes in structure. 




FIG. 19: Cross-section along the r-z plane of coherent vor- 
ticity of (2,0,1) ring mode. Top: (a) non-oscillated. Bottom: 
(b) oscillated. The vortex core shifts y + — 7 away from the 
wall. 



Thus, faster advection can be interpreted in two fash- 
ions. The first is in terms of the traveling wave. The 
shifting of the structures away from the wall into higher 
velocity mean flow causes these structures to travel 
faster, giving them less interaction time with the roll 
modes. Less interaction time with the roll modes means 
less energy transfer (less bursting), and due to their fast 
decaying nature, their lifetime is reduced. This reduced 
lifetime means they have less time to generate Reynolds 
stress, and therefore drag is reduced. The second inter- 
pretation is in terms of the classically observed hairpin 
and horseshoe vortices j 43 i 44 The pushing of the KL struc- 
tures away from the wall is equivalent to the vortices 
lifting and stretching away from the wall faster. This 
faster lifting and stretching process means that their life- 
time is shortened, again resulting is less time to generate 
Reynolds stress, and therefore drag reduction occurs. 



IV. CONCLUSIONS 

This work has shown, through a Karhunen-Loeve anal- 
ysis, four major consequences of spanwise wall oscillation 
on the turbulent pipe flow structures. They are: a shift- 
ing of rms velocities and Reynolds stress away from the 
wall; a reduction in the dimension of the chaotic attrac- 
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TABLE IV: Comparison of the measured location of the co- 
herent vortex core for two modes from each propagating sub- 
class in wall units away from the pipe wall . Each class 
shifts away from the wall, consistent with the shift in velocity 
rms and Reynolds stress. The (1,0,1) mode changes struc- 
ture significantly in the oscillated case as seen in Figure [TBI 
explaining the shift of its coherent vortex core towards the 
wall. 

Coherent vortex core location (wall units) 



Mode Non-oscillated Oscillated shift 



(1,2,1) 


41.8 


48.6 


6.8 


(1,5,1) 


28.7 


38.6 


9.9 


(2,2,1) 


45.6 


56.2 


11.0 


(3,2,1) 


65.9 


77.1 


11.2 


(1,1,1) 


47.2 


50.4 


3.2 


(2,1,1) 


53.9 


99.8 


45.9 


(1,0,1) 


31.9 


18.9 


-13.0 


(2,0,1) 


58.4 


65.8 


7.4 


(0,6,1) 


28.9 


39.0 


10.1 


(0,2,1) 


45.1 


52.9 


7.8 



tor describing the turbulence; a decrease in the energy 
in the propagating modes as a whole with an increase in 
modes that transfer energy to the outside of the log layer; 
and a shifting of the propagating structures away from 
the wall into higher speed flow resulting in faster advec- 
tion and shorter lifespans, providing less time to generate 
Reynolds stress and therefore reducing drag. 

The strength of the KL method is that it yields global 



detail and structure without conditional sampling. The 
ensemble was created out of evenly spaced flowfields in 
time, as opposed to conditional sampling of the nowfield 
with event detection such as bursts or sweeps, and the 
entire nowfield and time history was studied. Therefore, 
we argue that the overall mechanism of drag reduction 
through spanwise wall oscillation has been found. Al- 
though a result of drag reduction is the decorrelation of 
the low speed streaks and the streamwise vortices, as 
found by previous researchers, this is an incomplete de- 
scription of the dynamics. It is the lifting of the turbulent 
structures away from the wall by the Stokes flow induced 
by the spanwise wall oscillation that cause the reduction 
in the time and duration of Reynolds stress generating 
events, resulting in drag reduction. In addition, this dy- 
namical description encompasses other methods of drag 
reduction such as suction and blowing, active control, and 
ribblets, 22 establishing it as a contender for a universal 
theory of drag reduction. 
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